Ticks and Tick-borne Diseases
○ Elsevier BV
Preprints posted in the last 7 days, ranked by how well they match Ticks and Tick-borne Diseases's content profile, based on 12 papers previously published here. The average preprint has a 0.01% match score for this journal, so anything above that is already an above-average fit.
Joshi, K.; Susong, K. M.; Lim, A.; Liu, Y.; Brady, O. J.
Show abstract
Dengue is a mosquito-borne, viral disease of increasing public health significance. Currently, most public health interventions target the vector, with efficacy dependent on timing within the season. Whilst seasonal profiles have been characterised in some endemic settings a global assessment is lacking. Here, we develop and apply a proportion-based measure of dengue seasonality to reported case time series from 1990 to 2024 across 106 countries and territories, the largest assessment of this phenomenon to date. We identify regional differences in seasonality such that every month of the year saw cases peak in at least one country or territory. Latitude was identified as influencing seasonality, with cases peaking between March and April in the southern hemisphere and July and October in the northern hemisphere. Equatorial locations displayed flat seasonality, and amplitude increased with distance from the equator. K-means clustering identified three seasonal profile types: two with pronounced seasonal outbreaks (with distinct peak timing and shape) and one with flatter, more endemic transmission. Peak month timing covaried among locations within the same seasonality cluster, with phase differences meaning that information on shifts in peak timing may be available several months in advance in some settings, of potential significance for prediction and intervention planning. Beyond aiding public health planning, identification of seasonal clusters suggests that information on dynamics in one location could be leveraged to improve forecasting power in others with similar seasonal dynamics.
Ward, S.; Lawford, H.; Sartorius, B.; Mayfield, H.; Sam, F. A. L.; Sheridan, S.; Thomsen, R.; Viali, S.; Vaccher, S.; Robinson, L. J.; Angrisano, F.; Lau, C. L.
Show abstract
Background Serosurveillance can estimate the prevalence of antibodies (Ab) acquired through vaccination or current and/or past infection. Multi-pathogen serosurveillance that measures multiple Ab simultaneously, can enable identification of vulnerable sub-populations with immunity gaps to vaccine preventable diseases (VPD) and concurrent burden of neglected tropical disease (NTD), including those nearing elimination (lymphatic filariasis [LF], trachoma) and eradication (yaws). This study aimed to estimate seroprevalence and identify temporal trends of selected VPDs and NTDs in Samoa to inform targeted public health action. Methodology/Principal Findings Dried blood spots were collected from four repeated community-based surveys in eight primary sampling units (PSU) in Samoa in 2018, 2019, 2023 and 2024. Multiplex bead assays were used to detect Abs against antigens (Ag) for diphtheria, measles, rubella, tetanus, LF [Wb123 or Bm14], yaws [both Rp17 and TmpA; <14 years only], and trachoma [Pgp3; <14 years only]. Seroprevalence estimates were adjusted for sampling design and standardised for age and sex. Overall, 2,871 participants were included in this analysis. Seroprevalence of measles increased from 42% in 2018 to 95% in 2024, whereas diphtheria decreased from 79% in 2018 to 65% in 2024. Seroprevalence to yaws remained <1% for all years, whereas trachoma decreased from 21% to 7% (2018-2024). LF seroprevalence decreased between 2018 and 2024 for Bm14 (37% to 9%) and increased for Wb123 (10% to 22%). This study identified 15 (0.5%) individuals who were seronegative to all VPDs (7 in 2018; 8 in 2019); of these, five were seronegative to all VPDs and seropositive to at least one NTD. Conclusions/Significance Identification of sub-populations with concurrent seronegativity to VPDs and seropositivity to NTDs underscores the potential role of multi-pathogen serosurveillance in directing public health interventions to those at greatest risk. Examination of temporal patterns offer a valuable tool for measuring intervention impacts and progress towards elimination goals.
DA FONSECA, E. M.; Perry, K.; Barker, B.; Hirschi, M.; Hanson, K. E.; Walter, K. S.
Show abstract
Background Coccidioidomycosis is an emerging fungal disease across the arid Americas and a frequent cause of community-acquired pneumonia. Understanding where Coccidioides populations originate, how they move across space, and whether they are expanding is important for interpreting changing patterns of Valley fever and anticipating future infection risk. Methods We prospectively collected and whole-genome sequenced 186 Coccidioides-positive clinical isolates submitted to a national diagnostic laboratory, and included 126 previously sequenced genomes. We applied genomic clustering, time-calibrated phylogenetic reconstruction, ancestral area reconstruction, mating-type assignment, and demographic inference to identify major populations, infer dispersal patterns, assess evidence for recombination and clonality, and reconstruct historical population dynamics. Findings We analyzed 312 genomes (139 C. immitis; 173 C. posadasii) and identified three major genetic populations within each species. C. immitis included two California-centered populations and one Pacific Northwest population, whereas C. posadasii included two Arizona-centered populations and one Texas-centered population. The most recent common ancestor was estimated at approximately 127,000 years for C. immitis and 234,000 years for C. posadasii. Most populations were not fully monophyletic, consistent with retained ancestral variation and/or ongoing gene flow. Inferred dispersal was largely asymmetric, with most movement originating from California in C. immitis and from Arizona and Texas in C. posadasii. Most populations contained both mating types, but one C. immitis population and a Brazilian subgroup of C. posadasii were clonal. All populations showed recent demographic expansion. Interpretation The evolutionary history of Coccidioides is characterized by strong geographic structure, ongoing gene flow, and recent demographic expansion. These processes are likely to influence future patterns of Valley fever endemicity and supports the use of genomic surveillance to detect shifts in disease risk as environmental conditions change.
Bouhentala, O. W.; Kadir, M. Y.
Show abstract
Background. In 2026, the Democratic Republic of the Congo (DRC) experienced the largest recorded outbreak of Ebola disease caused by Bundibugyo virus, with epidemiologically linked importations and secondary transmission in Uganda. This study analysed the publicly reported trajectory and assessed the risk of introduction and onward transmission in North Africa and Europe. Methods. Public surveillance reports from the World Health Organization (WHO), European Centre for Disease Prevention and Control (ECDC), Africa CDC, national ministries of health, and peer-reviewed sources were synthesised through 15-17 July 2026. Headline counts and crude case-fatality ratios were restricted to laboratory-confirmed cases. Average notification rates were calculated from cumulative DRC counts. Exact Poisson intervals used the Garwood method, and the June-July rate ratio was estimated on the log scale. Risk was assessed across introduction likelihood, conditional onward-transmission likelihood, impact, and confidence. Results. By 15 July, the DRC had reported 2,124 confirmed cases and 828 deaths (crude confirmed-case fatality ratio, 39.0%) across 46 health zones in five provinces. Uganda had reported 20 confirmed cases and two confirmed deaths: 15 imported infections and five secondary cases, with no documented community transmission. DRC notifications averaged 47.4 per day during 1-15 July versus 35.9 per day during 2-29 June (rate ratio 1.32; counting-model 95% interval 1.20-1.46). WHO reported that more than 80% of new cases were detected outside known contact lists, while 119 confirmed healthcare-worker infections and 36 deaths had occurred. Introduction likelihood was assessed as very low to low for North Africa and very low for the general European population; delayed recognition in routine healthcare was the principal scenario for limited secondary transmission. Interpretation. Available indicators were inconsistent with effective control in eastern DRC at the data cut-off. Public reporting-date series cannot separate transmission from changing ascertainment, but they showed no sustained decline. Preparedness in North Africa and Europe should prioritise complete exposure histories, rapid isolation, validated diagnostics, protected clinical care, and contact management rather than reliance on border screening. Keywords: Bundibugyo virus; Ebola disease; outbreak surveillance; rapid risk assessment; importation; North Africa; Europe; Algeria; International Health Regulations.
Ma, Q.; Zhang, T.; Lin, D.; Zou, W.
Show abstract
Objectives: Although HIV incidence has declined in some settings, the overall global burden of sexually transmitted infections remains a major public health concern. In the context of the World Health Organization's call for people-centred STI prevention and care, identifying the shared geographic pattern of multiple STIs using data-driven analysis may help detect vulnerable areas and inform integrated prevention strategies. Methods: We analysed country-level incidence counts from the Global Burden of Disease 2023 study for 204 countries and territories over 1990-2023. A Bayesian shared-component spatiotemporal model was fitted, decomposing each disease's log-rate into a shared spatial component (scaled intrinsic conditional autoregressive prior), disease-specific spatial deviations, disease-specific first-order random walk temporal effects, and five socioeconomic covariates, with a negative binomial likelihood to accommodate overdispersion. The shared spatial score - the posterior mean of the shared spatial component - was used as a continuous index of STI co-occurrence burden. Posterior exceedance probabilities quantified directional stability. External validity was assessed via Spearman correlation with the Socio-demographic Index and generalised estimating equation regression of HIV/AIDS mortality on the shared score. Results: The shared spatial score exhibited marked geographic heterogeneity. The five highest-scoring countries were Eswatini (2.25), Lesotho (2.13), Malawi (1.90), Mozambique (1.89), and South Africa (1.85), all in southern Africa. Fifty-seven countries had high directional stability (posterior exceedance probability >0.95), concentrated in sub-Saharan Africa and the Caribbean. The score correlated negatively with SDI (Spearman rho = -0.619, p = 6.4 x 10^-23) and positively with HIV/AIDS mortality (incidence rate ratio = 14.64 per standard deviation, 95% CI: 11.90-18.01). Prior sensitivity analysis confirmed near-perfect ranking stability (rho >= 0.9999). Conclusions: STI co-occurrence is geographically concentrated, with the highest shared burden in sub-Saharan Africa and persistently elevated shared spatial signals also observed in parts of mainland Southeast Asia and the Caribbean. The shared spatial score provides a data-driven tool for prioritising integrated STI screening and prevention resources across countries.
Nyarko, E.; Antwi, P.; Amponsah, E. B.; Ofori-Boadu, L.; Oduro-Mensah, E.; Oliver-Commey, J. A.; Haruna, M.; Serwaa, C.; Dadzie, G.
Show abstract
Snakebite envenoming is a major neglected tropical disease disproportionately affecting rural populations in sub-Saharan Africa. In Ghana, evidence on the spatial and temporal distribution of risk remains limited, constraining targeted prevention and resource allocation. This study quantified district-level snakebite risk across Ghana, identified persistent hotspots and environmental drivers, and evaluated the relationship between snakebite burden and geographic access to treatment. Monthly district-level snakebite cases from Ghana's District Health Information Management System (2020 to 2025) were analyzed across all 261 districts using a Bayesian spatio-temporal model incorporating spatial effects, a temporal random effect, and a space-time interaction, fitted via Integrated Nested Laplace Approximation. Environmental covariates including rainfall, temperature, humidity, and NDVI quantified associations with risk. Relative risks, exceedance probabilities, Local Indicators of Spatial Association, and geographic accessibility identified priority districts. Snakebite risk showed strong spatial clustering and temporal variation. Persistent high risk districts were concentrated in Upper West (Daffiama Bussie Issa, Wa East, Wa West, Sissala East), Savannah (Bole, Gonja), North East (Mamprugu Moagduri), Western North (Bia East), Bono (Banda), Oti (Krachi Nchumuru), Western (Wassa East), and Eastern Region (Nsawam Adoagyiri, Fanteakwa North), though patterns evolved. Fanteakwa North emerged as the highest risk district nationally in 2025. Humidity and temperature were associated with increased risk, while rainfall and NDVI showed no significant effect. High risk districts often had poor treatment access, revealing inequities. This first nationwide Bayesian spatio temporal assessment provides an evidence base for surveillance, antivenom distribution, and interventions supporting WHO's 2030 snakebite reduction goals.
Bouhentala, O. W.; Kadir, M. Y.
Show abstract
Background. The 2026 Bundibugyo virus disease (BVD) epidemic in the Democratic Republic of the Congo (DRC) was declared on 15 May 2026 and determined a public health emergency of international concern on 17 May 2026. Public surveillance reporting consists of cumulative counts by report date; no line list with symptom-onset dates is available. Widely circulated characterisations - that this is the fastest-growing Ebola outbreak on record, that reported cases are doubling every 22 days, and that the case fatality ratio (CFR) is 37.5% - rest on these aggregates. We examined what the published data actually support. Methods. We assembled twelve published anchor points from 15 May to 13 July 2026 from WHO, WHO AFRO, NICD and the DRC National Institute of Public Health; one was recovered by back-calculation and checked against the directly reported subsequent total. We computed mean daily incidence between anchors and the within-interval death-to-case ratio. We reconstructed symptom-onset dates by Richardson-Lucy deconvolution with right-truncation correction under assumed onset-to-report delays with means of 5, 7 and 9 days, estimated the instantaneous reproduction number using the Cori method, and computed three CFR estimators: crude, resolved-case and outcome-delay-adjusted. Provincial CFRs used exact binomial intervals. Results. Confirmed cases plateaued at 40-52 per day for eighteen days, from 25 June to 13 July. Over the same period, the within-interval death-to-case ratio rose from 0.28 to 0.58. Reconstructed Rt was 1.28, with a 95% credible interval of 1.15-1.41, on 7 July, having fallen from approximately 2.9 in mid-May. Growth on the reconstructed onset curve corresponded to a doubling time of approximately 90 days, compared with 22 days computed from cumulative counts. CFR estimates were 37.5% for the crude estimator, 50.2% for the outcome-delay-adjusted estimator and 67.3% for the resolved-case estimator. Crude provincial CFR was 34.9% with a 95% confidence interval of 32.7-37.1 in Ituri, 58.2% with a 95% confidence interval of 50.7-65.5 in North Kivu, and 81.0% with a 95% confidence interval of 58.1-94.6 in newly affected provinces. Conclusions. A flat case count accompanied by a rising death-to-case ratio is difficult to reconcile with a transmission plateau and is consistent with saturated case detection. Reported case counts appear to have substantially decoupled from transmission and cannot presently distinguish control from detection failure. Doubling times computed from cumulative totals are artefacts. Test volume and positivity by health zone are the critical missing denominators.
Djaafara, B. A.; Elyazar, I. R.; Yosephine, P.; Surya, A.; Silalahi, F. S.; Handito, A.; Thohir, B.; Aryani, D.; Gunawan, D.; Nisa, A. K.; Prianto, E.; Samad, I.; Cook, A. R.; Huang, A. T.; Clapham, H. E.; Bhatt, S.; Mishra, S.
Show abstract
Estimating dengue force of infection (FOI) is essential for understanding transmission dynamics and targeting intervention programmes, yet surveillance data in endemic settings required for estimations are often incomplete, with varying formats. We developed a Bayesian hierarchical catalytic model that jointly fits age-stratified case data, aggregate case data, and seroprevalence surveys within a single framework, incorporating external covariates to improve parameter identifiability. Synthetic validation showed that covariates alone recovered accurate FOI point estimates even when most districts contributed only aggregate data, but did so with poorly calibrated uncertainty; anchoring the model with a single seroprevalence survey was necessary to bring credible interval coverage close to nominal. Applied to 128 districts across Java and Bali, Indonesia (2016-2024), the model revealed substantial spatial heterogeneity in FOI and reporting rates. Many districts in Java exceeded the WHO-suggested seroprevalence threshold for vaccine introduction, yet were classified as low-priority when using reported incidence as prioritisation criterion, particularly in areas with weak surveillance. Model-based seroprevalence estimation, integrating multiple data sources, offers a more consistent basis for identifying high-priority districts for vaccine introduction, and is less susceptible to surveillance bias than reported incidence.
Nimalrathna, S. U.; Harischandra, H.; Kimber, M.; Chandrasena, N.; De Silva, N.; Mallawarachchi, H.; De Silva, B. G. D. N. K.
Show abstract
The World Health Organization (WHO) validated Sri Lanka had eliminated lymphatic filariasis as a public health problem in 2016, the second country in Southeast Asia to attain this status. However, post-validation surveillance has identified sporadic cases of brugian filariasis. The reemergence of Brugia malayi infections in Sri Lanka warrants urgent investigations. Recent studies have shown that the parasite responsible for the reemergence is a novel zoonotic Brugia sp. maintained among dogs that is closely related but distinct to the human-infecting B. malayi species. The current study employed morphological and morphometric assessments, revealing that this novel zoonotic Brugia sp. is within the B. malayi morphological range. Molecular characterization of three genomic regions, the nuclear genomic region SLXI, the non-coding region HhaI, and the mitochondrial genomic region COXI confirmed it as a genetic variant more closely related to B. malayi than to B. pahangi. Phylogenetic analysis further indicated it as a distinct genomic variant, closely related to a B. malayi-like parasite reported from India. Notably, that same parasite was identified in infected humans, animals, and potential vector mosquitoes. This, together with the detection of both human and animal blood within the same brugian infective mosquitoes, and delineating the canine origin of the parasites in human infections, provides compelling evidence supporting zoonotic transmission of this parasite. To our knowledge, this is the first report demonstrating the presence of the same brugian parasite in humans, domestic animals, and potentially infective mosquitoes in Sri Lanka, supported by multi-genomic evidence. The recent identification of multiple potential mosquito vector species suggests that this parasite may have undergone adaptive changes, facilitating its ability to overcome the species barrier. These findings substantiate the long-held hypothesis of zoonotic transmission of the reemerged brugian parasite, highlighting significant implications for ongoing surveillance and control strategies.
Kim, S.; Mogasale, V. V.; Vesga, J. F.; Kang, H.; Skrip, L.; Jung, S.-m.; Islam, A.; Endo, A.; Edmunds, W. J.; Abbas, K.
Show abstract
Background Nipah virus (NiV) is a priority zoonotic pathogen causing high-fatality outbreaks. Early NiV outbreaks in Malaysia and Singapore had limited transmission beyond spillover events. However, since 2001, NiV outbreaks with person-to-person transmission have occurred in Bangladesh and India, driven by the NiV-Bangladesh genotype and NiV-India genotype. Our study aims to estimate the reproduction number, offspring dispersion, and serial interval governing NiV transmission in Bangladesh and India during 2001-2026. Methods We conducted a systematic review of NiV outbreak investigations in Bangladesh and India, searching PubMed, Embase, Web of Science, and grey literature through 28 February 2026. Case-level offspring counts from 27 eligible sources (323 cases across 67 outbreaks) were used as input to a hierarchical Bayesian negative binomial offspring distribution model. The serial interval was estimated by parametric distribution fitting to 137 transmission pairs. Country-stratified and sensitivity analyses were performed to evaluate the robustness of estimates. Results Pooling across 67 outbreaks, we estimated a median reproduction number of 0.46 (95% CrI: 0.28-0.73), an offspring dispersion parameter of 0.07 (0.05-0.10), and a serial interval of 13.3 days (95% CI: 12.8-13.8). Country-stratified median reproduction numbers were 0.48 (0.23-0.97) for India and 0.35 (0.19-0.59) for Bangladesh, and dispersion parameters were 0.04 (0.02-0.07) and 0.11 (0.06-0.18), respectively, indicating marked overdispersion in both settings. Conclusion NiV transmission is self-limiting on average and highly overdispersed, suggesting that a disproportionate share of onward transmission arises from a small number of cases. This epidemiological profile supports targeted containment measures, including contact tracing and quarantine, for effective NiV outbreak control.
Soubrier, H.; Seo, D.; Barks, P.; Meakin, S.; Mossoko, M.; Kitenge, R.; Dieberg, K.; Van Herp, M.; Mambula, C.; Flasche, S.; Camacho, A.; Coulborn, R.; Simons, E.; Ahuka-Mundeke, S.; Broban, A.
Show abstract
Background. The 2018-2020 Ebola virus disease outbreak in the Democratic Republic of the Congo (DRC) was the country's largest, and the second largest globally, amid armed conflict and community mistrust. Transmission heterogeneity (superspreading) is recognised in Ebola epidemics, but empirical estimates of its extent and determinants remain scarce for DRC outbreaks. We quantified transmission heterogeneity and its determinants during this outbreak. Methods. In this retrospective observational study, we reconstructed transmission chains for confirmed and probable cases (Aug 1, 2018, to June 25, 2020) using routinely collected Ministry of Health and Medecins Sans Frontieres surveillance data. We modelled the offspring distribution with a Bayesian negative binomial framework, correcting for incomplete contact tracing, to estimate the effective reproduction number (Reff), dispersion parameter (k), and proportion of cases responsible for 80% of transmission (prop80), overall, by subgroup, and over time. Individual-level determinants were assessed with a regression extension, adjusting for covariates. Findings. Among 3481 cases, 2008 transmission events linked 2402 (69%) individuals into 415 chains (median size 3, range 2-102). Overall Reff was 1.00 (95% CI 0.92-1.08) with k 0.29 (0.26-0.32); 17.8% of cases generated 80% of transmission. Overdispersion stayed stable despite fluctuating Reff. Non-isolation (IRR 1.79), death outside a treatment centre (IRR 4.34), and unfollowed contact status (IRR up to 4.48) predicted more secondary cases; vaccination cut transmission by about 60%. Interpretation. Epidemiological investigations linked 67.7% (2356/3481) of cases into 415 transmission chains (median size 3, range 2-102); linkage to a known infector fell to 10% during the November 2018-February 2019 period of peak insecurity. Transmission was heterogeneous overall, with dispersion parameter k of 0.29 (0.26-0.32), such that 17.8% (16.8-18.9) of cases generated 80% of onward transmission confirming superspreading as a stable, structural feature of Ebola dynamics. Critically, k remained stable throughout the outbreak, including during periods of elevated Reff, indicating that transmission surges reflected intensification of the same underlying process rather than new superspreading contexts, and that Reff alone is an insufficient summary of epidemic potential. Regression analyses identified predominantly modifiable determinants: cases not isolated in an Ebola treatment centre (IRR 1.79 [1.54-2.06]) or who died outside one (IRR 4.34 [3.47-5.31]) generated substantially more secondary cases, as did those registered as contacts but not followed up (IRR 3.01 [2.39-3.70]) or unregistered altogether (IRR 4.48 [3.67-5.38]) relative to actively followed-up contacts. Vaccination reduced onward transmission by 60-64% (IRR 0.36-0.40). These findings indicate that transmission was shaped less by gaps in epidemiological knowledge than by the operational reach of contact tracing, isolation, and vaccination delivery, particularly during periods of insecurity.
Pasaribu, A. P.; Nanine, I.; Ainur, F.; Jimanto, V.; Hutagalung, A. P.; Panggalo, L. V.; Devin, D.; Siregar, O. R.; Hasibuan, B. S.; Fahmi, F.; Trianty, L.; Coutrier, F. N.; Sasmono, R. T.; Satyagraha, A. W.
Show abstract
Red blood cell (RBC) disorders arose as an advantageous evolutionary response to malaria infections. In a heterozygous condition, such as in Southeast Asian Ovalocytosis (SAO), hosts are protected against severe malaria. In malaria-endemic regions, RBC disorders are presumed to be highly prevalent. Tanjung Leidong, a moderately endemic area in North Sumatra (API 1.13 in 2024), lacks comprehensive data on RBC disorder prevalence beyond G6PD deficiency. Therefore, this study aims to characterize the RBC disorders in this region as well as to characterize the anemia status in children living in Tanjung Leidong. Schoolers attending D. I. Panjaitan elementary to high school were recruited and screened for malaria by microscopy and G6PD deficiency using the STANDARD G6PD Assay. The DNA of the participants was also extracted to be genotyped for SAO, Hemoglobin E (HbE), and -thalassemia. Exclusively, G6PD-deficient DNA samples were genotyped further to determine variants. The proportion of G6PD deficiency, SAO, HbE, -thalassemia one-gene deletion, and two-gene deletion were 0.90%, 0.90%, 2.40%, 6.26%, and 0.30%, respectively. Anemia prevalence was approximately 14%, and RBC disorders were observed across children with normal to obese BMI. No malaria infections were detected by microscopy. The predominance of asymptomatic RBC disorders highlights that they are protective against malaria infection, although their protective role against malaria could not be directly assessed in this study. Both nutritional and genetic factors are found to contribute to anemia in this cohort. These findings underscore the importance of integrated screening strategies for RBC disorders and anemia in malaria-endemic settings.
Brochu, H. N.; Shi, Q.; Song, K.; Zhang, Q.; Munroe, J.; Harris, N. J.; Britt, N.; Zeng, Q.; Kapuria, K.; Chappell, J.; Norvell, B. M.; Peavy, L.; Williams, J. D.; Harris, A. B.; Chaitram, J.; Hutson, C. L.; Deng, J.; McGrath, D.; Boles, D.; Dale, S. E.; Gigante, C. M.; Iyer, L. K.
Show abstract
Background The 2022-2023 global mpox outbreak highlighted the critical need for robust genomic surveillance capabilities to track mpox virus (MPXV) evolution and transmission dynamics. Methods Building upon our established SARS-CoV-2 sequencing infrastructure, we implemented a Molecular Loop probe-based long-read sequencing approach using Pacific Biosciences Sequel II technology for comprehensive MPXV genomic surveillance across the United States (US). From August 2024 to June 2025, we generated 326 high-quality whole genome sequences from residual mpox-positive clinical specimens collected by Labcorp across all 10 US Department of Health and Human Services regions. Results Our analysis identified two samples containing clade Ib MPXV in January and June 2025 and captured shifting trends in clade IIb diversity, with 13 distinct lineages observed. We also identified multiple instances of large (~1.6-17.6kb) deletions proximal to the inverted terminal repeats in clade IIb genomes. APOBEC3 mutation analysis indicated substantial evidence of human-to-human transmission among both clades. Further, we observed significantly higher APOBEC3-associated SNPs per kilobase (P<0.001) in clade IIb genomic variable regions relative to their central conserved region. Our assay exhibited strong reproducibility across biological replicates from individual patients and accuracy was confirmed via parallel sequencing of select specimens by US Centers for Disease Control and Prevention (CDC) using metagenomic sequencing. We also demonstrated via custom simulation that our assay discriminates all known MPXV clades and lineages, including those we have not observed in the US. Conclusions Our integrated nationwide surveillance system facilitates real-time genomic tracking of outbreak evolution, with demonstrated capacity across SARS-CoV-2 and MPXV, positioning this platform for rapid deployment during future pathogen emergence.
Hussein, M. A.; Doshi, R.; He, L.; Reynolds, T.
Show abstract
Patients and caregivers seek informational and emotional support throughout medical care, especially when interpreting unfamiliar laboratory test results. Although resources such as patient portals and online health communities (OHCs) help address questions, gaps remain. The emergence of large language models (LLMs) offers the potential to be a complementary source of support to assist patients and caregivers in understanding and using their test results. The objective of our study is to empirically compare LLM responses to patients online questions containing their laboratory test results to responses written by peers in an OHC. We compared the 519 peer replies to 122 laboratory test-related posts from an OHC to 488 responses generated from four LLMs using mixed computational and qualitative methods. LLMs frequently provided clear explanations of medical terminology and structured interpretations of numeric results but were longer and less readable. Peers offered more personalized, context-specific emotional support. Overall, LLMs have the potential to complement peer responses in OHCs, but require greater emotional depth, reasoning transparency, and alignment with community norms.
Vijay, A.; Prabhune, A.; Srihari, V. R.; Rayampalli, A.
Show abstract
We present FootNet, a 453-image multi-view smartphone foot dataset for binary foot segmentation, with expertannotated masks across six anatomical views (dorsal, medial, and plantar, both left and right). We benchmark four segmentation models under a controlled protocol: U-Net with a MobileNetV2 encoder achieves the best performance (IoU 0.9268, Dice 0.9608, 95 % CI [0.9209, 0.9320]); DeepLabV3 with MobileNetV3-Large scores IoU 0.8984 (Dice 0.9449); UNet++ with MobileNetV2 scores IoU 0.8913 (Dice 0.9391); and SAM ViT-B with oracle boundingbox prompt scores IoU 0.9219 on the matched 191-image subset. Bonferroni-corrected Wilcoxon signed-rank tests (k = 6 comparisons) show U-Net significantly outperforms DeepLab (p < 0.001, r = 0.638) and SAM ViT-B with oracle boundingbox (p = 0.005, r = 0.202); UNet++ does not significantly differ from DeepLab (p = 0.062). Connected-component postprocessing yields negligible benefit (mean {triangleup}IoU = +0.0003, 12 of 453 images improved). The extended dataset is available upon request
Liu, J. B.; Chen, Y.-J.; Edelen, M. O.; Pusic, A. L.; Martin, N. E.; Zeng, C.
Show abstract
Purpose: Nonresponse to routinely collected patient-reported outcome measures (PROMs) threatens the representativeness of aggregated data. We characterized patient-, provider-, and clinic-level factors associated with PROMIS Global-10 nonresponse in routine radiation oncology care. Methods: In this retrospective cohort study, all adults seen at five Mass General Brigham radiation oncology clinics over one year were included. The primary outcome was patient-level nonresponse, defined as never completing the portal-administered Global-10 versus completing it at least once. Using iterative mixed-effects logistic regression, we modeled patient-, provider-, and clinic-level factors. Results: Among 12,214 patients, 71 providers, and five clinics, patient- and appointment-level response rates were 35.4% and 10.9%, with patient-level response ranging nearly fivefold across clinics (12.8% to 66.2%). In Model 1, male sex, lower education, not working, and recent surgery had higher odds of nonresponse, and longer time since diagnosis lower odds. After provider- and clinic-level factors were added, patient sex, education, and employment became nonsignificant, whereas recent surgery (adjusted odds ratio [aOR] 1.97) and longer time since diagnosis (aOR 0.46 for >12 months) persisted. A provider's historical collection rate was protective but attenuated at the clinic level. There, a later program launch (aOR 0.29) and higher historical collection rate (aOR 0.79) correlated with lower nonresponse, whereas academic versus community setting did not. Conclusions: Nonresponse to routinely collected PROMs is a multilevel phenomenon driven substantially by clinic-level implementation factors, not patient characteristics alone. Because response rate is only a proxy for representativeness, PROMs programs and PRO-based performance measures should prioritize representative collection over volume.
Weerasinghe, C.; Osowicki, J.; Simpson, J. A.; Crocker-Buque, T.; McCarthy, J.; Williams, E.; Price, D. J.
Show abstract
Controlled human infection models (CHIMs) are increasingly used in infectious disease research to study pathogen dynamics and evaluate interventions under controlled conditions. However, these studies are resource-intensive and involve ethical and safety constraints, making efficient study design critical. Dose-finding is a key early component in CHIMs, where the aim is to identify a challenge dose that achieves a target infection probability. Traditional rule-based designs are commonly used but can be inefficient, motivating the use of model-based adaptive approaches such as the Bayesian Continual Reassessment Method (CRM). Although CRM has been extensively studied and widely adopted in Phase I oncology trials for identifying the maximum tolerated dose of therapeutics, its application in CHIM settings remains limited, particularly when the endpoint of interest is infection. This tutorial provides step-by-step guidance for implementing a Bayesian CRM in dose-finding CHIMs, using an oropharyngeal Neisseria gonorrhoeae challenge as a motivating case study. The framework outlines key design components, including dose-grid specification, dose-response model, prior elicitation, Bayesian updating, decision rules, and stopping criteria, with particular emphasis on a clinically interpretable parameterisation. Trial operating characteristics are evaluated through simulation studies under multiple dose-response scenarios and prior-predictive analyses, and compared with a commonly used '3+3' type rule-based design. This work highlights the advantages of Bayesian model-based designs for dose-finding in CHIMs over classic rule-based designs and provides a structured, reproducible framework for implementing CRM, supporting their application in future CHIM studies.
Prosty, C.; Butler-Laporte, G.; Brophy, J.; Frenette, C.; Loo, V.; Coburn, B.; Hota, S.; Longtin, Y.; Kong, L.; Muller, M.; Steiner, T.; Valiquette, L.; Daneman, N.; Daley, P.; Nott, C.; MacFadden, D. R.; Kandel, C.; Chen, Y.; Perez- Patrigeon, S.; Lee, T. C.; McDonald, E.
Show abstract
Background and Aims The optimal treatment for first episodes and first recurrences of Clostridioides difficile infections (CDI) is unknown and there is emerging evidence for pulse and taper (P-T) regimens. Therefore, we sought to estimate the relative efficacy of treatment options. Methods MEDLINE and CENTRAL were searched from database inception to May 21, 2025 and unpublished conference abstracts were searched from recent infectious disease conferences. RCTs on the treatment of first episodes or first recurrences of CDI comparing fixed-dose or P-T regimens of fidaxomicin or vancomycin were included. The primary and secondary outcomes were 40- and 56-day CDI recurrence, respectively. A random-effects network meta-analysis on the risk ratio (RR) scale was conducted using a standard regimen (10-14 days) of vancomycin as the comparator. Treatments were ranked using the surface under the cumulative ranking curve (SUCRA). Results 8 RCTs were included comprising a total of 2181 patients. For 40-day recurrence, fidaxomicin P-T had the highest probability of ranking best (RR=0.10, 95%Confidence Interval [95%CI]=0.10-0.49, SUCRA=1.00), followed by vancomycin P-T (RR=0.49, 95%CI=0.32-0.76, SUCRA=0.61), fixed-dose fidaxomicin (RR=0.61, 95%CI=0.49-0.76, SUCRA=0.39), and, finally, fixed-dose of vancomycin (SUCRA=0.00). The treatments ranked in the same order for 56-day recurrence, though only 3 RCTs reported on this timepoint. Conclusion Vancomycin P-T, fidaxomicin P-T, and fixed-dose fidaxomicin were all superior to a fixed-dose vancomycin. Head-to-head comparative effectiveness RCTs are needed to quantify their relative effect sizes of and impact on long-term prevention of recurrent CDI.
Kamelian, K.; Pascall, D. J.; Cheng, M. T. K.; Meng, B.; Altaf, M.; Morse, R. M.; Aggio, J. B.; Egan, D. J. S.; Chen-Xu, M.; Trivioli, G.; Sutton, B.; Richter, A.; Gonzalez-Vazquez, L. D.; Cormie, C.; Kemp, S.; Yeadon, R.; Hyatt, B.; Wong, A.; Thesin Pelamkulangara, N.; Fraser, E.; McCarthy, B.; Novaes, F.; Stott, S.; Galvin, A.; Bellis, K. L.; De Angelis, D.; Harrison, E. M.; Martin, D.; Smith, R. M.; Gupta, R. K.
Show abstract
Background: Monoclonal antibodies have emerged as a prophylactic strategy to prevent symptomatic SARS-CoV-2 infection in immunocompromised individuals. However, the evolutionary and clinical implications of breakthrough infections under this regime remain unclear. Methods: A male in their 80s with a haematological/oncological diagnosis received a 2000 mg intravenous infusion of sotrovimab in March 2023 and was diagnosed with COVID-19 by RT-qPCR from a nasopharyngeal swab in August 2023. Weekly samples (n=24) were collected through February 2024 (171 days). All samples underwent whole-genome sequencing, with select mutations subjected to functional assessment. Findings: Sequencing identified the GE.1 lineage at all timepoints. An intra-host recombination event in ORF1ab (positions 8942-12458) was detected prior to 23 weeks post-detection, followed by a 14-fold increase in viral load (7.42e+06 to 1.00e+08 RNA copies/mL) and a marked shift in the viral population. E340D, a sotrovimab resistance mutation, was detected at low abundance (46%) within the first week post-infection, fluctuated over time, and was nearly fixed by week 15 (107 days) post-detection. We assessed five spike mutations - V36M, S98F, and V213G in the N-terminal domain, Y505P in the receptor-binding domain, and P681Q near the S1/S2 cleavage site - and additionally evaluated the impact of E340D. V36M conferred the highest infectivity across all cell lines, with the most significant effect in low-TMPRSS2 cells. While all mutations showed enhanced infectivity with the addition of E340D, the effect was most pronounced in mutations with lower baseline infectivity. The addition of E340D significantly decreased relative neutralizing titres for V36M, S98F, and V213G, enabling escape from neutralizing antibodies in XBB-responsive individuals, illustrating an enhanced phenotypic advantage. Patient neutralizing activity was absent pre-sotrovimab, and sotrovimab-induced neutralization was further compromised by selection of E340D. Interpretation: Sotrovimab pre-exposure prophylaxis in an immunocompromised patient did not prevent SARS-CoV-2 infection, and selected for resistant mutation E340D, with unexpected fitness consequences across non-receptor binding domain spike regions.
Moradi, E.; Dahnke, R.; Gaser, C.; Rikkonen, T.; Kroger, H.; Vaananen, S.; Solomon, A.; Sund, R.; Tohka, J.
Show abstract
Magnetic Resonance Imaging (MRI) derived brain age varies substantially between individuals, but it remains unclear whether early deviations from normal brain ageing precede future cognitive decline and whether they provide predictive value beyond conventional MRI measures. Here, we investigated whether MRI-derived brain age gap estimation (BrainAGE) identifies early structural brain ageing differences among cognitively normal individuals who later develop mild cognitive impairment (MCI) or dementia. We analysed longitudinal structural MRI data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) and replicated the main findings in the population-based Kuopio Osteoporosis Risk Factor and Prevention Study (OSTPRE). Individuals who later converted to MCI or dementia had higher BrainAGE values several years before diagnosis and, in ADNI, showed steeper longitudinal increases than stable individuals. Elevated BrainAGE values were also associated with increased risk of future conversion to MCI in cognitively healthy individuals and faster subsequent memory decline. Cross-sectional differences and the association between BrainAGE and risk of future conversion were replicated in OSTPRE. Importantly, adding BrainAGE to models including demographic, APOE4, cognitive, and MRI-derived measures consistently improved prediction of future cognitive outcomes, with the greatest benefit observed for individuals who converted after longer follow-up. These findings show that structural brain ageing begins to diverge years before the onset of MCI. BrainAGE captures this early divergence, providing complementary information beyond conventional structural MRI measures that may improve the early identification of cognitively normal individuals at increased risk of future cognitive decline when integrated with other biomarkers.